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The Revival of the Unified Dark Energy-Dark Matter Model ? 
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We consider the generalized Chaplygin gas (GCG) proposal for unification of dark energy and 
dark matter and show that it admits an unique decomposition into dark energy and dark matter 
components once phantom-like dark energy is excluded. Within this framework, we study structure 
formation and show that difficulties associated to unphysical oscillations or blow-up in the matter 
power spectrum can be circumvented. Furthermore, we show that the dominance of dark energy is 
related to the time when energy density fluctuations start deviating from the linear 5 ~ a behaviour. 

PACS numbers: 



I. INTRODUCTION 

The GCG model 0,0 is an interesting alternative to 
earlier proposals aiming to explain the observed accel- 
erated expansion of the Universe such as an uncanceled 
cosmological constant "s"] and quintessence P|, the latter 
being a variation of the idea that the cosmological term 
could evolve 0- 

In the GGG approach one considers an exotic back- 
ground fluid, described by the following equation of state 
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(1) 



where A and a are positive constants. The case a = 1 
corresponds to the Ghaplygin gas. In most phenomeno- 
logical analyzes the range < a < 1 has been consid- 
ered. Within the framework of Friedmann-Robertson- 
Walker cosmology, this equation of state leads, after in- 
serted into the relativistic energy conservation equation, 
to an energy density evolving as 
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(2) 



where a is the scale-factor of the Universe and B an in- 
tegration constant which should be positive for a well- 
defined Pch at all times. Hence, we see that at early times 
the energy density behaves as matter while at late times 
it behaves like a cosmological constant. This dual role 
is at the heart of the surprising properties of the GCG 
model. Moreover, this dependence with the scale-factor 
indicates that the GCG model can be interpreted as an 
entangled mixture of dark matter and dark energy. 
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The GCG model has been successfully confronted with 
various phenomenological tests: high precision Cosmic 
Microwave Background Radiation data Q, supernova 
data 0, and gravitational lensing In a recent work 
, it has been explicitly shown that regarding the latest 
supernova data [13, the GCG model is degenerate with 
a dark energy model involving a phantom-like equation 
of state (See also Ref for a detail study with differ- 
ent Supernova data sets). This excludes the necessity of 
invoking an unphysical fluid violating the crucial domi- 
nant energy condition for theoretical model building of 
our Universe which leads to the big rip singularity in 
future. The GCG, on the other hand, can mimic such 
an equation of state, but without any such pathology as 
asymptotically it approaches to a well-behaved de-Sitter 
universe. 

Despite all these pleasing features, the main concern 
with such an unified model is that it produces unphysical 
oscillations or even an exponential blow-up in the matter 
power spectrum at present p^ . This is expected from 
the behaviour of the sound velocity through the GCG 
fluid. Although, at early times, the GCG behaves like 
dark matter and its sound velocity is vanishing, as one 
approaches the present, the GCG starts behaving like a 
dark energy with a substantial negative pressure yielding 
a large sound velocity which, in turn, produces oscilla- 
tions or blow-up in the power spectrum. In any unified 
approach this is unavoidable unless one can successfully 
identify the dark matter and the dark energy components 
of the fluid. Naturally, these components are interact- 
ing as both are entangled within a single fluid. This is 
the main motivation of our present investigation. We 
show that the GCG is a unique mixture of interacting 
dark matter and a cosmological constant-type dark en- 
ergy, once one excludes the possibility of phantom-like 
dark energy. Due to the interaction between the compo- 
nents, there is a flow of energy from dark matter to dark 
energy. This energy transfer is vanishingly small until 
recent past, resulting in a negligible contribution from 
the cosmological constant at the time of gravitational 



collapse {zc — 10). This makes the model indistinguish- 
able from a CDM dominated Universe in the past. Just 
before the present (z ~ 2), the interaction starts to kick 
off producing a large energy transfer from dark matter to 
dark energy leading to its dominance at present. We have 
also shown that the epoch of this dark energy dominance 
is similar to that when dark matter perturbations start 
deviating from its linear behaviour. Moreover, in this ap- 
proach, the Newtonian equations for small scale pertur- 
bations for dark matter do not involve any fc-dependent 
term; hence, neither oscillations nor blow-up in the power 
spectrum are expected. We should mention that the de- 
caying dark matter model has been previously considered 
as an interesting possibility to solve, within the CDM 
model, the problem of over-production of dwarf galax- 
ies as well as the over-concentration of dark matter in 
halos 14]. Our results show that the GCG model is an 
interesting option for that scenario, such that the decay 
product is nothing but a cosmological constant. 
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FIG. 1: Qdm and Qa and Qt, as a function of redshift. We 
have assumed ildmo = -25, SIao ~ 0.7 and Q,i,o ~ 0.05 and 
a = 0.2. 



II. DECOMPOSITION OF THE GCG FLUID 



In Ref. |2I], it has been shown that the GCG La- 
grangian density has the form of a generalized Born- 
Infeld theory: 
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which clearly reproduces the Born-Infeld Lagrangian 
density for a = 1. The field 9 corresponds to the phase 
of a complex scalar field . 

Let us now discuss the decomposition of the GCG into 
components. Using Eqs. (^3), @ and introducing the 
redshift dependence, the pressure is given by 
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[yl + S(l + z)3(i+°)] i+« 
while the total energy density can be written as 
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where one has set the present value of the scale-factor, 
flo, to 1. 

Decomposing the energy density into a pressureless 
dark matter component, pdm, and a dark energy com- 
ponent, px with an equation of state wx, it follows that 
the equation of state parameter of the GCG can be writ- 
ten as 
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Thus, using Q), (jSJ and 0, one obtains for px 
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It is easy to see that, requiring that px > leads to the 
constraint wx < for early times {z ^ 1) and wx < —1 
for future (z ~ —1). Hence, one concludes that wx < —1 
for the entire history of the universe. The case wx < — 1 
corresponds to the so-called phantom-like dark energy, 
which violates the dominant-energy condition and leads 
to an ill defined sound velocity (see however 0). If one 
excludes this possibility, then the energy density can be 
split in an unique way: 
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and 
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from which one obtains the scaling behaviour of the en- 
ergy densities 
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One should note that this type decomposition based 
on the tachyon-like lagrangian has previously been con- 
sidered by Padmanabhan and Choudhury LLSj . Next we 
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express parameters A and B in terms of cosmological ob- 
servables. From Eqs. © and ((Tn|l . it follows that 



Pcha = PdmQ + PAo ^ (A + B) i+" , 



(12) 



where PchOj Pdmo ^nd pAo are the present values of Pchi 
Pm and Pa, respectively. Parameters A and B can then 
be written as a function of pcho 
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It is also useful to express A and B in terms of rJ^mOi 
r^Ao, the present values of the fractional energy densities 
^rfm(A) = Pm{A)/Pc whcrc pc is the critical energy density, 
Pc = SH^/SttG. Using the Friedmann equation 
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where p^o is the present baryon energy density, one ob- 
tains 



(15) 



Hence, with 



(22) 



Hence the evolution of dark energy and dark matter are 
linked so that energy is exchanged between these com- 
ponents (see Refs.0, 0| for earlier work on the inter- 
action between dark matter and dark energy). One can 
see from Figure 1, that until z ~ 2, there is practically 
no exchange of energy and the A term is approximatcUy 
zero. However, around z ~ 2, the interaction starts to 
kick off, resulting in an important growth of the A term 
at the expense of the dark matter energy density. Around 
z ~ 0.2, dark energy starts dominating the energy con- 
tent of Universe. Obviously, these redshift values depend 
on the a parameter and, in Figure 1, we have assumed 
a = 0.2. Nevertheless, the main conclusion is that in this 
unified model, the interaction between dark matter and 
dark energy is practically zero for almost the entire his- 
tory of the Universe making it indistinguishable from the 
CDM model. The energy transfer starts just in the recent 
past resulting in a significant energy transfer from dark 
matter to the A-like dark energy. In the next section we 
shall see that this epoch of energy transfer is similar to 
the one when dark matter perturbations start departing 
from its linear behaviour. 



III. 



STRUCTURE FORMATION 
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one can express the fractional energy densities Qdmi ^A 
and rjft as 
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where 
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Finally, given that Qdmo and f^Ao are order one quan- 
tities, one sees that at the time of nucleosynthesis, fl\ is 
negligibly small, making the model consistent with the 
nucleosynthesis process. 

Notice that there is an explicit interaction between 
dark matter and dark energy. This can be seen from 
the energy conservation equation, which in terms of the 
components can be written as 



In order to study structure formation, it is convenient 
to write the 0-0 component of Einstein's equation as 



= 8^G(pd,„ + pb) + A 
where A is given by 

A = SnGpA . 
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The energy conservation equation for the background 
fiuid is given by Eq. (|22|l . This is reminiscent of earlier 
work on varying A cosmology ^^, .18,, 19.] where the cosmo- 
logical term decays into matter particles. In our case, it is 
the opposite as a is always positive and hence the energy 
transfer is from dark matter to dark energy. This leads 
to the late time dominance of the latter and ultimately 
to the present accelerated expansion of the Universe. 

Let us now consider the issue of energy density per- 
turbations. We start by writing down the Newtonian 
equations for a pressureless fluid with background den- 
sity pdm , and density contrast Sdm with a source term due 
to the energy transfer from dark matter to the cosmolog- 
ical constant type dark energy. Assuming that both, the 
density contrast Sdm and peculiar velocity v are small, 
i.e. that Sdm << 1 and v << u, where u is the velocity 
of a fiuid element of volume, one can write the Euler, the 
continuity and the Poisson's equations in the co-moving 
frame as follows [T9| : 



FIG. 2: Sdm as function of scale factor. The solid, dotted, 
dashed and dash-dot lines correspond to a = 0, 0.2, 0.4, 0.6, 
respectively. We have assumed fldm ~ 0.25, fib = 0.05 and 
flA = 0.7. 
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FIG. 3: The growth factor m{y) as a function of scale factor 
a. The solid, dotted, dashed and dash-dot lines correspond 
to a = 0, 0.2, 0.4, 0.6, respectively. We have assumed fldm = 
0.25, fii, = 0.05 and = 0.7. 
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where $ is the gravitational potential, and is the source 
term in the continuity equation due to the energy transfer 
between dark matter and the cosmological constant type 
dark energy. The co-moving coordinate x is related to 
the proper coordinate r by r = ax. In our case, 
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Taking the divergence of Eq.(|2SJ) and using Eas. (PS|) 
and (|27|l . one obtains the small scale linear perturbation 
equation for the dark matter in the Newtonian limit: 
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Notice that, if = 0, i.e. without energy transfer, 
one recovers the standard equation for the dark matter 
perturbation in the ACDM case. One can also check that 
this occurs for a = 0. Moreover, one can easily see that, 
in the above equation, there is no scale dependent term 
to drive oscillations or blow up in the power spectrum. 

Let us turn to the evolution for the baryon pertur- 
bation in the Newtonian limit when the scales are well 



inside the horizon. Since we are considering the evolu- 
tionary period after decoupling, the baryons are no longer 
coupled to photons, there is no significant pressure due 
to Thompson scattering, and one can effectively consider 
baryon as a pressureless fluid like the dark matter fluid. 
Of course, the interaction between baryons and dark en- 
ergy is a relevant issue as it is related to the Equivalence 
Principle (see e.g. and references therein). In what 
follows we shall assume that the interaction of dark en- 
ergy with baryons vanishes. Thus, as baryons are uncou- 
pled, the strong interaction between the dark matter and 
dark matter violates the Equivalence Principle. Given 
that it is the parameter a that controls this interaction 
[a — Q means there is no interaction), it is a measure 
of the violation of the Equivalence Principle. One can 
also see from the behaviour of that this violation also 
starts rather late in the history of the Universe. In the 
Newtonian limit, the evolution of the baryon perturba- 
tion after decoupling for scales well inside the horizon is 
similar to the one described earlier for dark matter, but 
without the source term as there is no energy transfer to 
or from baryons. It is given by 
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where the third term in the l.h.s, we have dropped the 
contribution from baryons as it is negligible compared 
to the dark matter one. It is convenient to define for 
each component the linear growth function D(y) where 
y = log(a), 



D{y)5o 



(31) 



where 5o is the initial density contrast (assuming Gaus- 
sian distribution), as well as the growth exponent rn{y) = 
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FIG. 4: The bias 6 as a function of the scale factor, a. 
The solid, dotted, dashed and dash-dot lines correspond to 
Of = 0, 0.2, 0.4, 0.6, respectively. We have assumed Qdm = 
0.25, ^6 = 0.05 and Qa = 0.7. 

D'{y)/D{y). 

Asymptotically, as dark matter drives the evolution 
of the baryon perturbations, hence they grow with the 
same exponent m{y). However, their amplitudes may 
differ and their ratio corresponds to the so-called bias 
parameter, b = 5b/ 5dm- 

In what follows, we study the behaviour of 5drm "Tiiv) 
and b as function of the scale factor a. While solving 
the differential equations for the linear perturbation, the 
initial conditions are chosen such that at a = 10"'^, the 
standard linear solution 13 ~ a, is reached. In Figure 2, 
we have plotted the linear density perturbation for the 
dark matter, 6dmi as a function of a. One can see that 
while for a — (the ACDM case) the perturbation stops 
growing at late times, for models with a > 0, the pertur- 
bation starts departing from the linear behaviour around 
a ~ 0.8 (we have assumed for scale factor at present 
ao = 1) i.e. z ~ 0.25 which is similar to the epoch when 
the A term starts dominating (cf. Figure 1). In view of 
this behaviour, it is tempting to conjecture that, in our 
unified model, the interaction between dark matter and 
A-like dark energy is related with structure formation, so 
that sufficiently high density contrast {5dm >> 1) can re- 
sult a large energy transfer from the dark matter to the 
dark energy due to which acceleration of the Universe 
sets in. We should mention that this kind of scenario 
has been earlier discussed in Ref. I^J. Thus our study 
shows that there is a link between the structure forma- 
tion scenario and the dominance of dark energy which 
ultimately results in the acceleration of the Universe ex- 
pansion. This may give possible hint to the solution of 
so-called Cosmic Coincidence problem. 

The behaviour of the growth factor m(y) is also quite 
interesting. One can see from Figure 3, that between the 
present and z ~ 5, the growth factor is quite sensitive to 



the value of a. With a = 0.2, it increases up to 40% at 
present in relation to the value of the ACDM case. Notice 
that m{y) governs the growth of the velocity fluctuations 
in linear perturbation theory as the velocity divergence 
evolves as ^ = —Ham5dm] therefore, large deviations 
of the growth factor with changing a are detectable via 
precision measurements of large scale structure, through 
joint measurements of the redshift-space power spectrum 
anisotropy and bi-spectrum from z = to z ~ 2. 

We have also studied the behaviour of the bias param- 
eter in our model. In Figure 4, its evolution is shown. 
The plot suggests that the bias parameter also changes 
sharply in recent past with increasing a. This bias ex- 
tends to all small scales allowing for the Newtonian limit, 
being hence distinguishable from the hydrodynamical or 
nonlinear bias which takes place only for collapsed ob- 
jects. Thus, from the observation of large scale clustering 
one can distinguish the non-zero a case from the a = 
(ACDM) case. 

One can also see from Figure 2 that there is no sup- 
pression of 5dm at late times for any positive value of a, 
and thus one should not expect the corresponding sup- 
pression in the power spectrum normalization, erg, for the 
total matter distribution. This was one major problem in 
the previous GCG model approach which, as pointed out 
by Sandvik et al. cannot be solved even after the 

inclusion of baryons. In the new approach we propose 
here, one can successfully overcome this difficulty. 

Another interesting cosmological probe for our model 
comes from galaxy cluster M/L ratios. The most recent 
average value fi™ = 0.17 ± 0.05, has been determined 
by Bahcall and Comerford by observing 21 galaxy 
clusters around 2 ~ 1. The fact that nearby cluster data 
seem to prefer smaller values for Vim than the one ob- 
tained from WMAP data at z ^ 1100, can be regarded 
as a signal in favour of a decaying dark matter model like 
the GCG. 

Finally, one should expect a smaller ISW effect at early 
times, as the gravitational potential is practically con- 
stant up to z ~ 2 but, at late times, there will be a 
larger ISW effect due to the energy transfer from dark 
matter to dark energy resulting in a large variation of 
the gravitational potential. 



IV. DISCUSSION AND CONCLUSIONS 

In this work, we have considered a decomposition of 
the GCG in two interacting components. The first one 
can be regarded as dark matter since it is pressureless. 
The second one has an equation of state, px = t^xpx- 
It has been shown that ujx < — 1- Thus, once phantom- 
like behaviour is excluded our decomposition is unique. 
This shows that in a unified model like GCG one does 
not need anyhting more than a simple cosmological con- 
stant as dark energy which is consistent with the recent 
study by Jassal et. al 23] where it has been shown that 
a combination of WMAP data and observations of high 
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redshift supernovae, is fairly consistent with a cosmolog- 
ical constant like dark energy. 

It has been demonstrated that in this model the so- 
called dark energy dominance is related with the time 
when matter fluctuations become large {5dm > !)• Fur- 
thermore, we have shown that in what concerns structure 
formation, the linear regime {6dm ~ a) is valid fairly close 
to the present, meaning that at the time of structure for- 
mation begins, Zc ~ 10, the dark energy component was 
irrelevant and clustering occurred very much like in the 
CDM model. We have shown that the growth factor 
as well as the bias parameter have a noticeable depen- 
dence on the a parameter. We have implemented a model 



where there is a violation of the Equivalence Principle, as 
dark energy and baryons are not directly coupled. This 
may turn out to be a distinct observational signature of 
the present approach. 
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